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ABSTRACT 

We study the star formation history of normal spirals by using a large and homo- 
geneous data sample of local galaxies. For our analysis we utilise detailed models of 
chemical and spectrophotometric galactic evolution, calibrated on the Milky Way disc. 
We find that star formation efficiency is independent of galactic mass, while massive 
discs have, on average, lower gas fractions and are redder than their low mass counter- 
parts; put together, these findings convincingly suggest that massive spirals are older 
than low mass ones. We evaluate the effective ages of the galaxies of our sample and 
we find that massive spirals must be several Gyr older than low mass ones. We also 
show that these galaxies (having rotational velocities in the 80-400 km/s range) cannot 
have suffered extensive mass losses, i.e. they cannot have lost during their lifetime an 
amount of mass much larger than their current content of gas+stars. 
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1 INTRODUCTION 

The Star Formation Rate (SFR) is the most important and 
the less well understood ingredient in studies of galaxy evo- 
lution. Despite more than 30 years of observational and the- 
oretical work, the Schmidt law still remains popular among 
theoreticians and compatible with most available observa- 
tions (e.g. Kennicutt 1998a). 

It is well known that there are systematic trends in the 
SF history of the Hubble sequence: The ratio 9/ {9) of the 
current SFR \& to the past average one (\?) (integrated over 
the galaxy's age) increases as one goes from early to late type 
galaxies, albeit with a large dispersion within each morpho- 
logical type. The Hubble sequence seems to be determined 
by the caracteristic timescale for star formation, with early 
type galaxies forming their stars in shorter timescales than 
those of late types. However, this simple description of the 
Hubble sequence fails to answer two important, and proba- 
bly related, questions: i) what determines the caracteristic 
timescale for SF in galaxies? and ii) what is the role (if any 
at all) of a galaxy's mass? 

Gavazzi et al. (1996) found an anti-correlation between 
the SFR per unit total mass and galactic luminosity; part 
of this trend may reflect the aforementioned dependence of 
SFR on Hubble type, but it may also be that this trend 
is fundamentally related to the galactic mass. In a recent 
work Bell and de Jong (2000) use a large sample of spiral 
galaxies with resolved optical and near-infrared photometry 



and find that it is rather surface density that drives the star 
formation history of galaxies, while mass is a less important 
parameter. 

Photometric studies alone cannot lift the age- 
metallicity degeneracy, namely the fact that young and 
metal-rich stellar populations may be redder than old and 
metal-poor ones. Studies of the chemical aspects of galaxy 
evolution (i.e. gas fractions, star formation rates, metal 
abundances) can help to tackle the problem from a different 
angle, but they are also limited by the unknown history of 
gaseous flows from and to the system. A study combining el- 
ements from both photometric and chemical evolution offers 
the best chances to understand this complex situation. 

In this work we study the star formation history of 
spirals by using detailed models of chemical and spectro- 
photometric galactic evolution. The numerical code has been 
presented in detail elsewhere (Boissier and Prantzos 1999, 
hereafter Paper I) and is only briefly described in Sec. 2. 
It has been successfully applied to the modelling of global 
properties of spirals (Boissier and Prantzos 2000, hereafter 
Paper II), as well as to the corresponding abundance and 
photometric gradients (Prantzos and Boissier 2000, herefter 
Paper III). Based on the recent work of Boselli et al. (2000), 
we use here a large and homogeneous sample of data for 
normal spirals (presented in Sec. 3), ideally suited for the 
purpose of this work. Our results and the comparison to 
the observations are presented in Sec. 4. They suggest that 



© 2000 RAS 



S. Boissier, A. Boselli, N. Prantzos and G. Gavazzi 



galactic mass is the main driver of galactic evolution, al- 
though local surface density may also play a role. We argue 
that massive galaxies are, on average, older than less mas- 
sive ones, based on the fact that star formation efficiencies 
seem to be independent of galactic mass, while gas fractions 
are systematically low in massive spirals. Based on the ob- 
served star formation efficiencies and gas fractions, we de- 
rive in Sec. 5 corresponding galactic ages by using simple 
analytical models of galactic chemical evolution, taking into 
account gaseous flows. We find that low mass spirals must 
be several Gyr younger than massive ones. Also, we find that 
spirals with rotational velocities in the 80-400 km/s range 
must not have suffered extensive mass losses during their 
past history. Our conclusions are summarized in Sec. 6. 



2 MODELLING THE EVOLUTION OF DISC 
GALAXIES 

2.1 The model and the Milky Way 

In Boissier and Prantzos (1999, paper I), we presented a 
model for the chemical and spectrophotometric evolution of 
the Milky Way disc. We recall here the main ingredients of 
the model. 

The disc is simulated as an ensemble of concentric, in- 
dependently evolving rings, slowly built-up by infall of pri- 
mordial composition. For each ring we solve the classical 
equations of chemical evolution (e.g. Pagel 1997) without 
the assumption of Instantaneous Recycling Approximation. 
We use: stellar lifetimes from Schaller et al. (1992); the yields 
of Woosley and Weaver (1995) for massive stars and Renzini 
and Voli (1981) for intermediate mass stars; SNIa producing 
Fe with a rate given in Matteucci and Greggio (1986); and 
the Initial Mass Function of Kroupa et al. (1993) between 
0.1 M and 100 M . The adopted star formation rate (SFR) 
varies with gas surface density E 9 and radius R as: 



9(t,R) = aE g (t,R) 15 V(R) R~ 



(1) 



where V(R) is the rotational velocity, assumed to be 220 
km/s in the largest part of the Milky Way disc. This radial 
dependence of the SFR is suggested by the theory of star 
formation induced by density waves in spiral galaxies (e.g. 
Wyse and Silk 1989). The efficiency a of the SFR (Eq. 1) 
is fixed by the requirement that the observed local gas frac- 
tion, at R=8 kpc from the Galactic centre {o gas ~0.2), is 
reproduced at T=13.5 Gyr (our adopted value for the age 
of the local disc). 

The disc is built up by infall with a rate exponentially 
decreasing in time, and a characteristic time-scale r in f(R) 
increasing with radius (as to mimic the inside-out formation 
of the disc). The relation between T in f(R) and S(7?) (total 
surface density) is shown on Fig. 1 and is a posteriori justi- 
fied, since it is crucial in shaping the various radial profiles 
(of gas, SFR, abundances, etc) of the disc, which compare 
favourably to observations. 

The spectro-photometric evolution is followed in a self- 
consistent way, with the metallicity dependent stellar tracks 
from the Geneva group (Schaller et al. 1992, Charbonnel 
et al. 1996) and stellar spectra from Lejeune et al. (1997). 
Dust absorption is included according to the prescriptions 
of Guiderdoni et al. (1998) and assuming a "sandwich" con- 
figuration for the stars and dust layers. 



It turns out that the number of observables explained 
by the model is much larger than the number of free pa- 
rameters. In particular the model reproduces present day 
"global" properties (amounts of gas, stars, SFR, and su- 
pernova rates), as well as the current disc luminosities in 
various wavelength bands and the corresponding radial pro- 
files of gas, stars, SFR and metal abundances; moreover, 
the adopted inside-out star forming scheme leads to a scale- 
length of ~4 kpc in the B-band and ~2.6 kpc in the K-band, 
in agreement with observations (see paper I). 



2.2 Scaling relations 

For a simplified extension of the model to the the case of 
other disc galaxies we adopt the "scaling properties" de- 
rived by Mo, Mao and White (1998, hereafter MMW98) in 
the framework of the Cold Dark Matter (CDM) scenario 
for galaxy formation. For the details, the reader should re- 
fer to MMW98 and Boissier and Prantzos (2000, paper II). 
Discs form inside non baryonic dark matter haloes of var- 
ious masses and concentration factors. Assuming constant 
disc to halo mass ratios (here taken to be 771^=0.05), discs 
are characterized by two parameters: Vc, the circular veloc- 
ity measuring the mass of the disc, and A, the spin parameter 
measuring its angular momentum. A disc is described by its 
scale- length Rd and its central surface density Eo, which can 
be related to the ones of our Galaxy (designated by MW) 
by: 
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The total mass of the disc Md is proportional to Vc, but in- 
dependent of A. The A distribution deduced from numerical 
simulations (see MMW98 and references therein) presents 
a maximum at A ~ 0.04 — 0.05 and extends from ~0.01 
to ~0.20. Since the value of Xmw is not far from the peak 
value, we explored here values of A in the range 1/3 Xmw to 
3 Xmw- Assuming that A and Vc are independent, we con- 
structed models with velocities in the observed range 80 to 
360 km/s. We calculated the velocity profile resulting from 
the disc plus an isothermal dark halo, and we used Equ. 1 to 
calculate self-consistently the SFR (with the same coefficient 
a, "calibrated" on the Milky Way). 

At this point it should be noticed that our scaling re- 
lations (2 and 3) are based on the assumption that angular 
momentum is conserved during the evolution of the disk. 
Numerical hydrodynamical simulations of galaxy formation 
(e.g. Navarro and Steinmetz 1999 and references therein) do 
not support this idea; indeed, it is found that as the baryonic 
halo gas cools down and collapses to the disk it loses most of 
its angular momentum. In those conditions, the final config- 
uration of the disk cannot be related in a simple way to the 
initial A of the halo. However, such simulations lead , in gen- 
eral, to disk sizes much smaller than observed. It is possible 
that star formation and feedback are not properly described 
in those simulations. For instance, Sommer-Larsen, Gelato 
and Vedel (1999) found that delaying the cooling of the gas 
reduces the loss of angular momentum. Since the situation is 
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not clear yet, we make here the assumption of disk evolution 
at constant angular momentum, a posteriori justified by the 
fact that the resulting disk sizes are in agreement with ob- 
servations (see Paper II). For a more detailed discussion of 
these issues, see also Cole et al. (2000). 



2.3 Infall timescales 

The infall rate in our models is exponentially dependent on 
time, with a time-scale r in f that depends on the surface 
mass density and is calibrated on the Milky Way disc (see 
Fig. fil); as explained in paper II, we found that a depen- 
dence on the total mass of the galaxy is necessary in order 
to reproduce the observations of present-day discs. We con- 
sider that all the galaxies started forming their stars 13.5 
Gyr ago and that the formation of their exponential discs 
(characterized by Eo and Rd) was completed at the present 
epoch. Notice that the value of T=13.5 Gyr plays no essen- 
tial role in the results of this work (values in the 10-15 Gyr 
range would lead to quite similar conclusions); these results 
depend mainly on the infall timescales. The more massive 
galaxies are characterized by shorter formation time-scales, 
while less massive galaxies are formed on longer ones (Fig. 
hi). This assumption turned out to be a crucial ingredient 
of our models, allowing to reproduce an impressive amount 
of observed properties of spirals that depend on mass (or 
Vc)'- colours, gas fractions, abundances and integrated spec- 
tra. Those observables are thoroughly presented in paper II, 
while some of them are revisited in this work. In Prantzos 
and Boissier (2000, paper III in the series), we show that 
this model also reproduces fairly well the observed colour 
and abundance gradients in spirals. 

We notice that Bell and de Jong (2000) suggested re- 
cently that the observed colour gradients of disc galaxies cor- 
relate very well with the local surface density, in the sense 
that inner and denser regions are older. Our assumption 
about the infall timescale agrees, at least qualitatively, with 
their findings. 

We stress that infall timescales are inputs to the model, 
and adjusted as to reproduce observations. Star formation 
timescales are outputs of the model (resulting from the 
adopted prescriptions for infall and SFR) and are presented 
in Sec. 4.1. 



3 COMPARISON TO OBSERVATIONS 

3.1 The observational sample 

The sample of galaxies analysed in this work, which has been 
extracted from the large multifrequency database of nearby 
galaxies of Gavazzi and Boselli, is extensively described in 
Boselli et al. (2000) . Here we give just a brief description of 
the sample selection criteria: we refer the reader to Boselli 
et al. (2000) for the detailed references on the data and on 
their analysis. 

Galaxies analysed in this work are taken from the 
Zwicky catalogue (CGCG, Zwicky et al. 1961-1968) (m pg < 
15.7). They are either late-type (type>S0a) members of 3 
nearby (recession velocity cz < 8000 km s _1 ) clusters (Can- 
cer, A1367, Coma), or located in the relatively low-density 
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Figure 1. Infall time-scale r in f as a function of the local sur- 
face density S. The dependence on £ is calibrated on the Milky 
Way (solid curve, V _ c=220 km/s). The dependence on Vc is cho- 
sen as to reproduce the properties of the homogeneous sample of 
spiral galaxies presented in Sec. 3. When Ti n f < 0, the infall is 
exponentially increasing with a time-scale It-^jI. 



regions of the Coma-A1367 supercluster (ll h 30 m <RA< 
13 h 30 m ; 18°<dec<32°) as defined in Gavazzi et al. (1999a). 
To extend the present study to lower luminosities, we include 
in the sample the late-type Virgo cluster galaxies brighter 
than nipg < 14.0 listed in the Virgo Cluster Catalogue as 



cluster members (VCC, Binggeli et al. 1985). Furthermore 
VCC galaxies with 14.0 < m pg < 16.0 included in the "ISO" 
subsample described in Boselli et al. (1997a) and CGCG 
galaxies in the region 12 h <RA< 13 h ; 0°<dec<18° but out- 
side the VCC, are considered. 

To avoid systematic environmental effects we consider 
the subsample of late-type galaxies whose HI deficiency (de- 
fined as d = \og(H I / (H I)) , the ratio of the HI mass to the 
average HI mass in isolated objects of similar morphologi- 
cal type and linear size, see Haynes & Giovanelli 1984) is 
d >0.3, typical of unperturbed, isolated galaxies. The final 
combined sample comprises 233, mainly "normal" galaxies. 

We assume a distance of 17 Mpc for the members (and 
possible members) of Virgo cluster A, 22 Mpc for Virgo clus- 
ter B, 32 Mpc for objects in the M and W clouds (see Gavazzi 
et al. 1999b). Members of the Cancer, Coma and A1367 clus- 
ters are assumed at the distance of 62.6, 86.6 and 92 Mpc 
respectively. Isolated galaxies in the Coma supercluster are 
assumed at their redshift distance adopting H = 75 km s _1 
Mpc" 1 . 

For the 233 optically selected galaxies, data are avail- 
able in several bands as follows: 100% have HI (1420 MHz) 
and 99% H band (1.65 /im) data, while a much coarser cov- 
erage exists in the UV (2000 A)(29 %), CO (115 GHz)(38%) 
and Ha (6563 A) (65%), as shown in Table 1 and 2 of Boselli 
et al. (2000) . The distribution of the sample galaxies in the 
different morphological classes is given in Table 3 of Boselli 
et al. (2000). 
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As previously discussed, the present sample is optically 
selected and thus can be biased against low surface bright- 
ness galaxies; the inclusion of the Virgo cluster should in 
principle favor the presence of some low surface brightness 
galaxies, not easily detectable at higher distances. Being vol- 
ume limited, the sample is not biased toward bright, giant 
spirals, but it includes also dwarfs and compact sources. Its 
completeness at different wavelengths makes this a unique 
sample suitable for statistical analysis. 



3.2 Data analysis 

Ha and UV fluxes, corrected for extinction (and [Nil] con- 
tamination) as described in Boselli et al. (2000), are used to 
estimate star formation rates through population synthesis 
models given in that work; a power-law IMF of slope -2.5 
with a lower and upper mass cutoff of 0.1 Mq and 80 Mq, 
respectively, is adopted. Its high mass part is quite similar 
to the one of the IMF of Kroupa et al. (1993), adopted in 
our models. Given the uncertainty in the UV and Ha flux 
determination, in the extinction correction and in the trans- 
formation of the corrected fluxes into SFRs via population 
synthesis models, we estimate an uncertainty of a factor of 
~3 in the determination of the SFR. 

The total gas content of the target galaxies, HI+H2, 
has been determined from HI and CO measurements. CO 
(at 2.6mm) fluxes have been transformed into H2 masses as- 
suming a standard CO to H2 conversion factor of 1.0 10 20 



mol cm (K km 



(Digel et al. 1996). For galaxies 



with no CO measurement, we assume that the molecular 
hydrogen content is 10% of the HI, as estimated from iso- 
lated spiral galaxies by Boselli et al. (1997b). The total gas 
mass has been corrected for He contribution by 30%. HI 
fluxes are transformed into neutral hydrogen masses with 
an uncertainty of ~ 10%. The average error on CO fluxes is 
~ 20%; the error on the H2 content, however, is significantly 
larger (and difficult to quantify) due to the poorly known 
CO to H2 conversion factor (see Boselli et al. 1997b). 

Galaxy colours have been determined from broad band 
near-IR and optical photometry. Near-IR (H band) images 
are available for 230 of the 233 sample galaxies, while B 
images or aperture photometry for 214 objects. H and B 
magnitudes have been corrected for extinction as in Gavazzi 
& Boselli (1996). No correction has been applied to galaxies 
of type later than Scd. The estimated error on B and H 
magnitudes is 15 %. 

Rotational velocities have been determined from the 
HI line width at 21cm, and corrected for inclination as in 
Gavazzi (1987). To avoid large systematic errors, we esti- 
mate rotational velocities only for galaxies with inclinations 
> 30 deg and with the 21cm line width accurately deter- 
mined (double or single horned profile with high signal-to- 
noise). The uncertainty on the determination of the rota- 
tional velocity is ~ 15 km s _1 . 

In summary, we have a homogeneous sample of 233 nor- 
mal (non-perturbed) disc galaxies. For 96 of them, we eval- 
uated all the quantities of interest in this work: Blue magni- 
tude Mb, total mass Mt, total gas mass M g , star formation 
rate ^ and rotational velocity Vc- We shall see below how 
our models fit those properties and what kind of inferences 
can be made on the star formation history of those galaxies. 
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Figure 2. Colour index B — H vs rotational velocity Vc- Sym- 
bols are observations presented in Sec. 3.1. A clear correlation is 
obtained, with massive discs being redder than low mass ones. 
The curves are the results of our models and are labeled by the 
value of the ratio A/Amvc, where Xmw is the spin parameter for 
the Milky Way disc. The thin curve corresponds to very low val- 
ues of A, producing very compact galaxies that look like bulges 
or ellipticals rather than spirals. 



3.3 Mass-driven colours and gas fractions 

In Fig. 2 we compare our model results to the colours B — H 
of our galaxy sample. A clear correlation is obtained between 
B — H and Vc, with the more massive discs being, on aver- 
age, redder than their lower mass counterparts. Our models 
(solid curves) reproduce naturally the observed correlation, 
since by construction the more massive discs form their stars 
earlier (Fig. 1). The corresponding timescales for star forma- 
tion will be discussed in Sec. 4.1. The observed dispersion in 
the lower panel of Fig. 2 can be accounted for by the range of 
A values in our models, but only partially. Discs with lower A 
are more compact, have higher central surface densities and 
evolve earlier than those with larger A. Obviously, A values 
larger than 3 Xmw (the largest value used here), would lead 
to even smaller B — H values than the models displayed on 
Fig. 2, possibly accounting for the rest of the scatter. 

We notice that our models do take extinction by dust 
into account (with the prescriptions presented in Sec. 2.1), 
and that extinction contributes somewhat (by ~0.5 mag) 
in redenning the most massive discs in our models, which 
have relatively large amounts of metals. However, we insist 
on the fact that it is age, not extinction, which is mainly 
responsible for the trends of our models on Fig. 2. 

In Fig. 3 we present two important pieces of data, on 
which the main argument of this paper is based. In the upper 
panel, we show the ratio M g j Lh (mass of gas to H-band 
luminosity) vs rotational velocity Vc- Since Lh reflects the 
mass of the stellar population (Gavazzi et al. 1996), this 
ratio is a measure of the gas fraction of the system. Despite 
the large scatter, it is clear that a correlation does exist, the 
more massive galaxies being, in general, gas poor. Our model 
results can fairly well describe the observations, accounting 
for both the slope and the dispersion of the M g / Lh vs Vc 
relation. 

In the lower panel of Fig. 3 we plot the ratio \&/Lh (cur- 
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Figure 3. Upper panel: Ratio M g /Lu of gas mass to H— band lu- 
minosity vs rotational velocity Vc', it is a measure of gas fraction. 
Lower panel: Ratio ^ / Ljj of star formation rate to //—band lu- 
minosity vs. rotational velocity Vc', it is a measure of the present 
SFR compared to the past average one. Massive discs have smaller 
gas fractions and were more active in the past than their lower 
mass counterparts. In both panels, our models (parametrised with 
the spin parameter X/Xmw, where \mw 1s the Milky Way value 
of A) are in fair agreement with the data. 



rent SFR to H— band luminosity) vs. Vc- As discussed in 
several places (e.g. Kennicutt 1998a), this ratio is a measure 
of the parameter 6 = ^//{^/), the ratio of the current SFR 
to the past average one. Massive galaxies display smaller 
\P /Lh, and thus lower b values, than their low mass coun- 
terparts. This means that they formed their stars at much 
higher rates in the past. Again, our results compare fairly 
well with the data, concerning both the slope and the scat- 
ter of the correlation. For a given Vc, large A discs have 
larger gas fractions today and were less active in the past 
than their lower A counterparts. As we shall see in Sec. 4.3, 
and discuss in Sec. 5, when the two panels of Fig. 3 are com- 
bined, they suggest quite convincingly that low mass spirals 
are on average younger than massive ones. 

At this point, we notice that Gavazzi and Scodeggio 
(1996) already advanced a similar hypothesis to explain the 
observed colors of galaxies as a function of their H-band 
luminosity (a mesure of their dynamical mass according to 
them). Gavazzi et al. (1996) came to similar conclusions and 
suggested that mass is the main parameter of galaxy evolu- 
tion, on the basis of multiwavelength observations concern- 
ing a variety of disc properties (colours, gas content, star 
formation rate, radius, surface brightness). 
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Figure 4. Time-scale of an exponential fit to the star formation 
rate history of our model galaxies; T E ~ xp (Gyr^ 1 ) is plotted as a 
function of disc rotational velocity Vc and parametrised with the 
spin parameter X/\mw- F° r l° w values of Vc, Texp is negative, 
i.e. the star formation rate is increasing with time. i"^, xp is larger 
for massive galaxies, which form their stars earlier. 



4 STAR FORMATION HISTORIES IN DISCS 

4.1 Formation time-scales 

The Hubble sequence of galaxies is usually interpreted in 
terms of different star formation timescales (e.g. Kennicutt 
1998a), although such an interpretation leaves unclear the 
role of the galaxy's mass (Prantzos 2000). In order to put 
our results in that context, we performed an exponential fit 
to the star formation histories of our model galaxies (ex- 
cluding the first 2 Gyr, where such a fit turns out to be 
inadequate). The resulting timescales are shown in Fig. m 
Notice that in some models the SFR is continuously increas- 
ing with time, resulting in negative caracteristic timescales. 
For that reason we present Tg XP , which has the advantage 
of varying continuously when going from a SFR increasing 
in time to one decreasing with time. For an exponential star 
formation rate, ^ oc exp~ t ' TEXP , t^ xp is equal to — ^/^ 
(where ^ = dty/dt) and can be considered as the normalized 
rate of change of the SFR. Obviously, the larger is t^ xp , the 
earlier the galaxy forms its stars. 

As can be seen on Fig. 4, massive galaxies in our models 
form their stars on shorter timescales than their lower mass 
counterparts. This general trend is somewhat modulated by 
the spin parameter A: discs with smaller A (i.e. more com- 
pact) have larger T^ X p than discs with larger A of similar 
rotational velocity. The most massive discs of our simula- 
tions (Vc=360 km/s) have decreasing SFR with texp ~4-5 
Gyr. Discs with Vc ~200 km/s have long timescales, of the 
order of ~10 Gyr, i.e. essentially constant SFR. Finally, low 
mass discs ( Vc ~100 km/s) have SFR increasing in time, 
with caracteristic timescales r^i 



-0.2 to -0.5 Gyr" 



4.2 The evolution of the star formation efficiency 

What is the reason for the vastly different SF timescales 
obtained in our models as a function of Vc? Is the overall 
SF efficiency e = \& jM g directly affected by the mass of the 
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Figure 5. Evolution of the Star Formation Efficiency e (ty/M g ) 
as a function of time for three families of models labeled by their 
A values (1/3, 1, and 3 times Xmw)- Low A values correspond 
to more compact galaxies, while values larger than 3 \mw lead 
to Low Surface Brightness Galaxies. For each A, three curves are 
shown corresponding to Vc=80, 220 and 360 km/s, respectively. 
The SF efficiency e depends more on A than on Vc- At T=13.5 
Gyr, the present-day observations (Sec. 4.3) are shown; an artifi- 
cial small age spread is introduced. Our grid of models at T=13.5 
Gyr covers well the observed e values. 



galaxy? An inspection of Eq. 1 suggests that this cannot be 
the reason. Indeed, at the caracteristic radius of the disc R d , 
the local efficiency is: 



e(R d 



^{R d )/T, g {R d 



a £r V R-, 



A 



-2 ^0.5 ^ 



since R d oc A V from Eq. 2 and T,(R d ) oc Eo oc X~ 2 V from 
Eq. 3. In other terms, the SF efficiency at R d varies very little 
with Vc and depends much more on A than on Vc- Since the 
value of any intensive quantity (like SF efficiency) at R d is 
typical of the whole disc, it is expected that the global SF 
efficiency also depends little on Vc- In order to show this 
quantitatively, we plot in Fig. 5 the evolution of the global 
SF efficiency e of our models as a function of time. We show 
the results for 3 values of A (1/3, I and 3 times Xmw) an d 
three values of Vc (80, 220 and 360 km/s, respectively). 
We compare our results at T=13.5 Gyr to estimates of e = 
$ /M g in our sample of nearby spirals (presented in Sec. 
4.3) . The following points should be noted concerning our 
models: 

i) The efficiency e depends very little on Vc, especially 
during the last half of the history of the galaxy. 

ii) e is mainly determined by A: "compact" galaxies have 
higher efficiency because of their smaller size (^ oc 1/R, for 
a given Vc) and larger gas surface density at the caracteristic 
radius R d . At T=13.5 Gyr, the variation of our e values due 
to A can fully account for the dispersion in e measured in 
nearby spirals. 

iii) The SF efficiency in compact galaxies (low A) 
presents a peak at early times and then decreases. This hap- 
pens because the star formation migrates to outer regions 
(due to the inside-out disc formation scheme) where the lo- 
cal SF efficiency is lower (because of lower surface densities 
and of the I/R factor). In more extended galaxies (larger A) 



e does not present such a decrease because local properties 
vary little with radius. 

iv) For the lowest disc velocities Vc, e may increase 
with time. This is due to the adopted form of infall: the gas 
surface density increases considerably when the gas arrives 
finally in the disc (which may take a very long time in the 
case of the largest A values and lowest Vc values, see Fig. 

4). 

The main point of this section is that in our models 
the SF efficiency e does not depend directly on the mass 
of the galaxy. The range of e span by our models during 
galactic evolution is due to A and tinf, not to Vc', the most 
important of the two parameters is tinf- Of course, tinf is 
adjusted to Vc, so that the observations of Fig. 2 and 3 (and 
many others, presented in Paper II) are reproduced. This is 
a crucial ingredient for the success of our models, and we 
discuss its implications in Sec. 5. 

4.3 Gas fraction and star formation efficiency 

The analysis of Sec. 4.1 and Sec. 4.2 lead to an important 
conclusion: the success of our models is to be interpreted 
in terms of mass-dependent SF timescales; but this is not 
due to any explicit dependence of SF efficiency on galaxy 
mass. Indeed, the SF efficiency of our models is virtually 
independent of mass during most of galactic history, and in 
particular at the present time. Is this supported by observa- 
tions? 

In Fig. 6 we display our data of Fig. 3, this time in a 
more "physical" presentation, appropriate for a quantitative 
discussion (see Sec. 5). The gas fractions a gas (upper panels) 
and SF efficiency e (lower panels) are presented as a function 
of rotational velocity Vc (left panels), H-magnitude (middle 
panels) and B — H colour (right panels) . 

The gas fraction a ga3 = M g /Mr (where the total mass 
Mt = M, + M g is the mass of stars + gas) is obtained 
by converting the Lh luminosity to star mass M, through 
the mass to light ratio M„JLh obtained in our models. This 
value is in the 0.3-0.6 range in solar units (M,/Lh ~0.35 for 
V c =80 km/s, ~0.47 for V c 220 km/s and ~0.56 for Vc=360 
km/s). Despite the uncertainties, a clear trend is present 
in the upper panel: the more massive, luminous and red a 
galaxy is, the smaller is its gas fraction. Notice that a sim- 
ilar trend is obtained by McGaugh and de Blok (1997) and 
Bell and de Jong (2000). Our model results are in excellent 
agreement with the data, although we have some difficulty 
in reproducing blue and gas poor discs (with B — H <2.5 
and a gas <0.4). 

The SF efficiency of our observed galaxies is obtained 
by simply dividing the SFR ^ with the gas mass M g . We 
notice that the uncertainties in deriving e are rather large: 
those concerning the SFR ^ are quite large (a factor of ~3), 
whereas those of the gas mass at least of ±20 %. As a result, 
the observationally derived scatter of e, as appears on Figs. 
5 and 6, is certainly larger that the real one. It is clear, how- 
ever, that the SF efficiency does not seem to be correlated 
with the mass or the colours of spirals. The scatter in the 
observed values of e (a factor of ~10), should be compared 
to the range of ~60 span by galaxy mass (since the mass 
of the disc m d oc Vc) or the range of five magnitudes in 
H-luminosity. Our model SF efficiencies e at T=13.5 Gyr, 
also displayed in Fig. 6 (lower panel, solid curves), show no 
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Figure 6. Gas fraction (cr 9a3 = M s /Mt, upper panels) and global Star Formation Efficiency (e = ty/M g , lower panels) vs circular 
velocity Vq (left panels), H-band magnitude (middle panels) and colour index B-H (right panels). Dots represent our data (Sec. 4.3). 
The results of our models (Sec. 4.2) are shown by curves, parametrised by the the ratio A/Amw, where \mw is the Milky Way spin 
parameter. The gas fraction decreases, on average with galaxy's mass, luminosity and colour index, while the global SF efficiency e does 
not depend on those parameters. Our models account fairly well for these data. We note that e mot j e i is a function of A and is greater for 
compact galaxies. 



dependence on mass or colour. They also reproduce the ob- 
served dispersion, the more compact discs (smaller A) being 
the more efficient in turning their gas into stars. 

The upper and lower panels of Fig. 6, when combined, 
point to an important conclusion: since the SF efficiency is 
independent of galactic mass (or luminosity), the fact that 
low mass galaxies have larger gas fractions today may only 
be due to their smaller ages. This is the most straightfor- 
ward interpretation, independent of any theoretical consid- 
erations (except for the implicit assumption that the SF 
efficiency has remained ^constant during the galaxy's his- 
tory). This conclusion is corroborated by an independent 
observable, namely that low mass galaxies are, in general, 
bluer than more massive ones. Notice that the latter ob- 
servable concerns also elliptical galaxies, but the well known 
problem of the age-metallicity degeneracy does not allow to 
conclude in that case. In the case of spirals, the situation 
is even worse in principle, because colours may be affected 
by the presence of dust (presumably more abundant in mas- 
sive spirals). Because of this complication, the observed gas 
fraction (smaller in large spirals) is not sufficient in itself 
to lift the age-metallicity degeneracy. However, when com- 
bined with the fact that the SF efficiency is independent of 
galactic mass (as argued here), the degeneracy is lift, and 
the aforementionned conclusion is naturally obtained. 



Before turning to a more detailed discussion of our find- 
ings, we would like to point out that similar results are ob- 
tained in other recent works. For instance, McGaugh and 
de Blok (1997) find a clear trend between the gas fraction 
and the B— magnitude of their galaxy sample, that they de- 
scribe by the relation a 9as =0.12 (Mb+23). This relationship 
is shown in Fig. 7 (upper panel, dashed line), along with the 
McGaugh and de Blok (1997) data for normal spirals and 
our data for normal spirals. It can be seen that there is a 
very good agreement between the two data sets. We notice 
that our models fit the observations of McGaugh and de 
Blok (1997) fairly well (see Fig. 13 of Boissier and Prantzos 
2000). 

In the lower panel of Fig. 7 we show the SF efficiency, 
both for our data set and the one of Kennicutt (1998b). No- 
tice that Kennicutt (1998b) gives the average SFR surface 
density (\I f ) and average gas surface density (S 9as ) of nor- 
mal spirals, i.e. the integrated quantities are divided by the 
disc surface aerea (within the optical radius) . Obviously, the 
ratio (iff) / (Egas) gives the overall SF efficiency e (since the 
disc aerea cancels out). As can be seen in Fig. 7, Kennicutt's 
values of the SF efficiency are slightly larger than ours (by a 
factor of ~ 2), and have a smaller dispersion in the adopted 
logarithmic scale. Taking into account the various uncer- 
tainties in estimating the SFR from the data (see Sec. 3), 
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Figure 7. Gas fraction (a ga s = M 9 /Mt, upper panels) and 
global Star Formation Efficiency (e = \I//M 9 , lower panels) vs 
B— band magnitude. Our data are shown by filled symbols in both 
panels. They are in agreement with data from other surveys (Mc- 
Gaugh and de Blok 1997, squares in the upper panel; and Ken- 
nicutt 1998b, triangles in the lower panel). The same trends as 
in Fig. 6 are also obtained here (i.e. gas fraction decreasing with 
galaxy's luminosity and Star Formation efficiency independent of 
luminosity). The dashed curve in the upper panel is a fit to the 
McGaugh and DeBlok (1997) data: o- sas =0.12(Af s +23). 



such a discrepancy between Kennicutt's results and ours is 
not unexpected. But the important point is that Kennicutt's 
values are also independent of the galaxy B— luminosity and, 
by virtue of the Tully-Fisher relation, on the galaxy's mass. 
The anti-correlation between gas fraction and luminos- 
ity in spirals was noticed by several authors (e.g. Gavazzi 
et al. 1996, McGaugh and de Blok 1997, Bell and de Jong 
2000, Boselli et al. 2000, etc.). Based on two independent 
samples, we showed here that the SF efficiency of spirals is 
independent of their mass. The two findings combined, point 
to small discs being younger, on average, than massive ones. 
Our model, presented in Sec. 3.2, nicely explains these fea- 
tures (and several other, presented in papers II and III). 
However, one may argue that the complex interdependence 
between the adopted infall and SFR prescriptions makes a 
straightforward interpretation difficult; he/she may also ar- 
gue that other types of models could also account for the 
observations and give a different interpreration for the same 
data (e.g. by invoking outflows). For that reason, we discuss 
in the next section this issue on the basis of simple analyt- 
ical models of galactic chemical evolution and we make a 
very rough evaluation of the "ages" of the galaxies in our 



sample. Our purpose is not to derive the exact ages of the 
galaxies, but rather to have an order of magnitude estimate 
and, in particular, to check whether there is any trend of the 
derived ages with galactic mass. 



5 GALACTIC AGES 

In the framework of simple models of galactic chemical evo- 
lution adopting the Instantaneous Recycling Appproxima- 
tion (IRA), one may obtain analytical solution for various 
quantities. In particular, provided that the Star Formation 
Rate ^ is proportional to the gas mass M g , i.e. 



# = eM 



(5) 



one may obtain a relationship between the gas fraction a gas 
and time T (assuming that the SF efficiency e is constant 
in time). The form of this relationship depends on further 
assumptions about the evolution of the system, i.e. on the 
possibility of allowing for gas flows in or outside the "box" 
(e.g. Pagel 1997). 

Assuming that the galaxies of our sample have evolved 
as simple, homogeneous, "boxes" , we consider three possibil- 
ities: a "closed box" (all the gas is present from the very be- 
ginning), an "infall" model (where gas mass is continuously 
added to the system) and an "outflow" model (with gas con- 
tinuously leaving the system). In the cases of gaseous flows, 
further assumptions about the corresponding flow rates are 
required in order to obtain analytical solutions. More specif- 
ically: 
a) Closed Box: In that case, we have 



= exp[-(l - R)eT] 



(6) 



where the return fraction R accounts for the gas returned 
by stars to the interstellar medium; for the IMF of Kroupa 
et al. (1993) adopted here we have R ~0.32. 
b) Infall: An analytical solution is easily obtained if it is 
assumed that the infall rate just balances the gas deple- 
tion due to star formation. As we have shown in Paper II, 
with detailed numerical models reproducing a large body of 
observational data, this situation describes rather well the 
largest period in the lifetime of most spiral galaxies. In that 
case, we have: 



= [l + e(l -R)T] 



(7) 



c) Outflow: Analytical solutions may be obtained by as- 
suming that the outflow rate is proportional to the SFR 
font — 7 1 3/. In that case, we have: 

R 



1 



(R - l)exp[-e(i? - 1 - i)T\ - 7 



(8) 



We shall assume here that the outflow rate is equal to the 
SFR, i.e. 7=1 (since, for higher outflow rates, we obtain 
ridiculously low galactic ages). 

For each of the three scenarios, we shall consider five 
combinations between the observed gas fraction and SF ef- 
ficiency, in order to derive galactic age T through Eqs. 6, 7 
and 8. For the gas fraction vs. rotational velocity, we adopt 
the three curves of the upper panel of Fig. 8, correspond- 
ing roughly to the mean trend and the upper and lower 
bounds of the observations, respectively. For each of those 
three curves, we derive the corresponding galactic ages by 
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Figure 8. Upper panel: Observed gas fraction vs. galaxy's ro- 
tational velocity; the three curves represent, schematically, the 
mean trend ( thick curve) and the upper and lower bounds of the 
observed gas fractions in Fig. 6. Given the Star Formation ef- 
ficiency e, the shaded aerea inside the upper and lower curves 
may be used to derive approximate ages for the galaxies, in the 
framework of simple analytical models for galactic chemical evolu- 
tion, using IRA (Instantaneous Recycling Approximation). This 
is done in the three next panels, for a Closed Box, a model with 
Infall (such as the gas mass remains constant in time) and a model 
with Outflow (with an outflow rate equal to the Star Formation 
Rate). In each of those panels, the thick curve corresponds to the 
mean gas fraction of the upper panel and a SF efficiency e=0.3 
Gyr -1 (the mean value for the efficiencies in Fig. 6). The shaded 
aerea is obtained by keeping the same SF efficiency, but consid- 
ering the upper and lower bounds for the observed gas fraction in 
the upper panel (lower gas fractions lead to larger ages for a given 
rotational velocity). The dotted curves are obtained with the low 
gas fractions combined to a high SF efficiency (e=l Gyr - 1 , a rea- 
sonable upper bound for the observed e values in Fig. 6). Finally, 
the dashed curves are obtained with the high gas fractions com- 
bined to a low SF efficiency (e=0.1 Gyr - 1 , a reasonable lower 
bound for the observed e values in Fig. 6). 



the upper panel of Fig. 8) , combined to a high SF efficiency 
e=l Gyr -1 (a reasonable upper bound for the observations 
of Fig. 6 and 7) ; in that way we obtain the dashed curves in 
Fig. 8. 

An inspection of the results, plotted in the three lower 
panels of Fig. 8, shows that the derived "effective" galactic 
age is a monotonic function of rotational velocity Vc in all 
cases. The absolute ages depend, of course, on the adopted 
gas fractions, SF efficiencies and assumptions about gaseous 
flows. For a given gas fraction and SF efficiency, the infall 
model leads to the largest ages; since the gas is constantly 
replenished, it takes more time to attain a given gas frac- 
tion than in a closed box (starting with <r gas = l). Also, the 
outflow model produces the lowest ages; since part of the 
gas is constantly removed, it takes less time to reach a given 
gas fraction than in a closed box. The ages derived in the 
framework of the outflow model are, in general, too low: they 
are lower than 6 Gyr for all galaxies with Vc <200 km/s. 
It transpires that galaxies in this velocity range should not 
have suffered extensive mass losses, i.e. they must not have 
lost during their lifetime an amount of gas as important as 
their stellar content (since we adopted f ou t — $ here). 

The mean age values in the case of the infall model are 
4-5 Gyr for the small discs (in the Vc ~100 km/s range) and 
~10 Gyr for Vc ~220 km/s. For discs with Vc >300 km/s, 
extremely large ages (>15 Gyr) are found; however, the ap- 
proximation of evolution at constant gas mass is certainly 
not valid in that case (see Fig. 4 in Boissier and Prant- 
zos 2000). We notice that Bell and de Jong (2000), on the 
basis of a different sample and with a completely different 
method (based on a photometric estimate of the ages) find 
an "effective age" of ~10 Gyr for the more luminous galaxies 
(Mk ~-26) and ~6 Gyr for the less luminous ones (Mk ~- 
20); the galaxies of our sample span also this luminosity 
range, and its interesting to see that similar ages are found 
with completely independent methods. 

The results obtained in this section confirm the sugges- 
tion of Sec. 3.2: low mass discs are, on average, younger than 
massive ones. The only way to reverse this trend is by as- 
suming that the SF efficiency is strongly correlated with disc 
mass. Despite the large scatter in the observational data, it 
is clear that such a correlation does not exist, at least at the 
present time; and, since the observed galaxies span a large 
range in masses and metallicities, it does not seem plausible 
that such a correlation ever existed in the past. In the case 
of outflow models, another possibility would be to consider 
that the more massive discs suffered more important mass 
losses (leading to low gas fractions without having to invoke 
large ages). However, such a hypothesis is incompatible to 
the fact that massive discs have deeper potential wells and 
are less prone to outflows than low mass ones. 



using a SF efficiency e=0.3 Gyr -1 (the mean value of the 
observations in Fig. 6 and 7); in that way, we obtain the 
shaded regions in Fig. 8, with the thick curve corresponding 
to the mean trend, the largest ages to the lowest gas frac- 
tions and vice versa. We perform two more calculations. For 
the first, we adopt the high gas fractions (uppermost curve 
in the upper panel of Fig. 8), combined to a low SF efficiency 
e=0.1 Gyr -1 (a reasonable lower bound for the observations 
of Fig. 6 and 7) ; in that way we obtain the dotted curves in 
Fig. 8. Finally, we adopt the low gas fractions (lower curve in 



6 SUMMARY 

In this work we investigate the properties of the star for- 
mation efficiency of spiral galaxies and study the implica- 
tions for their evolution. We use a large homogeneous sam- 
ple of disc galaxies, for which we measure gaseous mass M g 
(HI), rotational velocity Vc, star formation rates ^ and lu- 
minosities in the B and H bands. We are then able to de- 
rive the corresponding gas fractions a gas and SF efficiencies 
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e = W /M g as a function of Vc, H— luminosity and B — H 
colour index. We find that the gas fraction is correlated to 
Vc, H— luminosity and B — H , in the sense that more mas- 
sive, luminous and redder discs have smaller gas fractions. 
Previous work by McGaugh and de Blok (1997) reached sim- 
ilar conclusions. The main finding of this work is that the 
SF efficiency does not correlate with any of the galaxy prop- 
erties; despite a rather large dispersion (within a factor of 
^10), the observed e is independent of Vc, Lh or B — H . 

We interpret our data in the framework of detailed mod- 
els of galactic chemical and spectrophotometric evolution, 
utilising metallicity dependent stellar lifetimes, yields, tracks 
and spectra. These models are calibrated on the Milky Way 
disc (paper I) and use radially dependent star formation 
rates, which reproduce observed gradients (paper Iff). They 
are extended to other spirals in the framework of Cold Dark 
Matter scenarios for galaxy formation, and are described by 
two parameters: rotational velocity Vc and spin parameter 
A. As in our previous work (paper II), we find good agree- 
ment with the observations, provided a crucial assumption 
is made: massive discs are formed earlier than less massive 
ones. With this assumption our models reproduce the ob- 
served trends of gas fractions and SF efficiencies vs. Vc, 
while variations due to A account for the observed disper- 
sion in both cases. 

It is important to notice that the dependence of age on 
galactic mass that we find is not due to any mass-dependent 
SF efficiency, only to disc formation timescales; in our mod- 
els this is achieved by varying the infall timescales. Both 
observations and models suggest that the SF efficiency is in- 
dependent of galaxy properties. Since the observed galaxies 
cover a wide range of masses, colours and metallicities, there 
is no reason to suppose that the SF efficiency was different 
in the past. The adopted SFR prescription in our models 
also results in a very slowly varying SF efficiency with time. 

When the observed relations of gas fraction vs Vc and 
SF efficiency vs Vc are considered in combination, they 
convincingly suggest that low mass discs are, on average, 
younger than more massive ones; this conclusion is indepen- 
dent of any model and the only assumption is that the SF ef- 
ficiency is constant in time (a quite plausible assumption, as 
argued above) . In the framework of simple analytical mod- 
els of galactic chemical evolution, we evaluate the "effective 
ages" for the galaxies of our sample, using the observation- 
ally derived values of a gas and e. We find that even models 
with modest outflows (with ejected masses equal to the stel- 
lar ones) lead to ridiculously low values for the galaxy ages; 
our conclusion is that galaxies in the range Vc ^80-400 km/s 
have not suffered extensive mass losses. Closed box models 
and infall models lead to more plausible values for the effec- 
tive ages. In particular, infall models lead to ages of ~4-5 
Gyr for discs of Vc ~100 km/s and ~8-10 Gyr for Vc ~200 
km/s. These "chemically derived" ages are in fair agreement 
with those derived on the basis of our more sophisticated 
numerical models, which fit a much larger body of observa- 
tional data for low redshift spirals. Most importantly, they 
are also in fair agreement with the "photometric" ages de- 
rived in a completely independent way and with a different 
sample by Bell and de Jong (2000). 

In summary, our data, taken at face value, suggest that 
the bulk of stars in more massive discs are older than in less 
massive ones. This is supported by our detailed numerical 



models of galactic chemical and photometric evolution, but 
also by recent, independent, analysis (Boselli et al. 2000). 
We notice that Bell and deJong (2000) conclude that it is 
local surface density that mainly drives the star formation 
history, while mass plays a less important role. Our study 
suggests that mass is the main factor, while local surface 
density plays only a minor role (through the spin parameter 
A: lower A values lead to higher local surface densities for a 
given rotational velocity Vc)- 

We notice that this picture is hardly compatible with 
the currently popular "paradigm" of hierarchical galaxy for- 
mation, which holds that large discs are formed by merging 
of small units at relatively late epochs. If this were the case, 
massive discs should have large SF efficiencies, in order to 
have their gas fractions reduced to lower levels than their 
less massive counterparts. However, such an enhanced SF 
efficiency is not supported by observations. 
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